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1 Introduction 

Now is a golden age of cosmology and astrophysics. Many abstractive notions such as black 
holes, dark matter (DM), dark energy etc. have become step by step more scientifically 
feasible and widely accepted subjects. According to the WMAP [1] the non baryonic dark 
matter, which is called the cold dark matter (CDM), must exist and contain approximately 
22% of all energy density of the universe. The characteristic of the CDM is thought to 
be massive and rarely interact with ordinary matter. There is yet little astrophysical data 
which bear on the CDM. However, there are a few proposals to explain the DM in the 
context of particle physics [2]. The most popular particles in this class are the sterile 
neutrino, the axion, the lightest supersymmetric particle and etc. 

The DM does not exist within the SM. It has to be realized beyond the SM at the 
electroweak scale or above, so that newly introduced particles in those models are poten- 
tially good candidates for the DM. In most supersymmetric models, there is a conserved 
multiplicative quantum number R- parity, which implies that the lightest superpartner is 
stable and can be a DM candidate. This kind of model can not only explain the origi- 
nal DM but also represent the greatest expectations in particle physics at TeV scale to be 
probed by the LHC. However, there is yet no experimental evidence to support the models. 
The other way to extend the SM is the enlargement of the gauge symmetry group from 
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SU(3)c® SU(2)i(gi U(l)y to larger groups. In particular, there exists a simple extension 
of the SM gauge group to SU(3)c'8) SU(3)L(g) U(l)x, the so called 3-3-1 models [3, 4]. 

Depending on the electric charge of particle at the bottom of the lepton triplet, the 
3-3-1 models are classified into two main versions: the minimal model [3] with the lepton 
triplet (z^, I, and the version with right-handed (RH) neutrinos [4], where the RH neu- 
trinos place at the bottom of the triplet: {u, I, i^'^)^. In the 3-3-1 model with right-handed 
neutrinos, the scalar sector requires three Higgs triplets. It is interesting to note that two 
Higgs triplets of this model have the same U(l)x charges with two neutral components at 
their top and bottom. In the model under consideration, the new charge X is connected 
with the electric charge operator through a relation 

Q = Ts-^Ts + X. (1.1) 

Assigning these neutral component vacuum expectation values (VEVs) we can reduce the 
number of Higgs triplets to two. Therefore, we have a resulting 3-3-1 model with two 
Higgs triplets [5, 6]. As a consequence, the dynamical symmetry breaking also affects the 
lepton number. Hence it follows that the lepton number is also broken spontaneously at a 
high scale of energy. Note that the mentioned model contains a very important advantage, 
namely, there is no new parameter, but it contains very simple Higgs sector; therefore, 
the significant number of free parameters is reduced. To mark the minimal content of the 
Higgs sector, this version that includes right-handed neutrinos (RHzy) is going to be called 
the economical 3-3-1 model. 

We would like here to emphasize that by choosing different electric charge operators, 
we can get a few different versions of 3-3-1 model such as the minimal 3-3-1 model, the 
3-3-1 model with right handed neutrino, the economical 3-3-1 model and etc. However, all 
those models have the same motivations such as 

1. The family number must be a multiplicative of three. 

2. It could explain why the value sin^ 6w < j is observed. 

3. It can solve the strong CP problem. 

4. It is the simplest model that includes bileptons of both types: scalar and vectors 
ones. 

5. The model has several sources of CP violation. 

Besides those motivations, the 3-3-1 models can contain a candidate for the dark matter. As 
an example, the 3-3-1 model with right handed neutrino exhibites that the scalar bilepton 
is a candidate for the dark matter as shown in Ref. [7], in which the authors provided a 
new quantum number to forbid the interaction of bilepton with the SM particles. However, 
in the economical 3-3-1 model, which we are considering now, the model naturally contains 
a candidate for the dark matter even without introducing a new quantum number or a 
discrete symmetry. 
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The paper is organized as follows: Section 2 is devoted to a brief review of the model. 
The Higgs sector is considered in the effective approximation w ^ v,u m section 3. Here, 
we also find out stable Higgs - a candidate for the dark matter. The relic abundance as 
well as its dependence/indepedence of parameters is figured out in section 4. We study 
direct and indirect searches for dark matter in section 5. Conclusions are given in the last 
one - section 6. 

2 A brief review of the model 

The particle content in this model, which is anomaly free, is given as follows: 

■0aL = (z^aL,^aL, (l^ai?)'')^ ~ (3, -1/3), /aij~(l,-l), a = 1,2,3, 

QlL = {uiL, diL, Uif ~ (3, 1/3) , 

QaL = (doL, -■UaL,-DaL)^ ~ (3*,0), Z^aij ~ (1 , - 1/3) , a = 2, 3, 

UaR ^ [1,2/2,) , daij- (1,-1/3), C/r~ (1,2/3), (2.1) 

where the values in the parentheses denote quantum numbers based on the (SU(3)l, U(l)x) 
symmetry. Unlike the usual 3-3-1 model with right-handed neutrinos, where the third 
family of quarks should be discriminating [8], in this model under consideration the first 
family has to be different from the two others [9] . The electric charges of the exotic quarks 
U and Da are the same as of the usual quarks, i.e., qu = 2/?,, qo^ = —1/3. 

The spontaneous symmetry breaking in this model is obtained by two stages: 

SU(3)l ® U(l)x ^ SU(2)i ® U(l)y ^ U(1)q. (2.2) 

The first stage is achieved by a Higgs scalar triplet with a VEV given by 

X = {xi,X2,xlV ~ (3,-1/3) , (X) = ^ {u,0,uf. (2.3) 

The last stage is achieved by another Higgs scalar triplet needed with the VEV as follows 

0= (<A^0^,</.^)^~(3,2/3), {<P) = ^{0,v,0f . (2.4) 

The Yukawa interactions which induce masses for the fermions can be written in the 
most general form: 

= ^L'NC + -^LNV, (2.5) 

in which, each part is defined by 

ClNC = h^QlLXUR + h^pQaLX*DpR 

pmn 

+hiQlL^daR + hl^QaL<P*UaR + H.C., (2.6) 
-^^LNV = S^QlLX'^aR + SaaQaLX*daR 

+s^QiL^DaR + s'^QaL(p*UR + H.c., (2.7) 
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Table 1. Nonzero lepton number L of the model particles. 



Field 








X? 


X2 


03 




DaL,R. 


L 


1 


1 


-1 


2 


2 


-2 


-2 


2 



Table 2. ;B and £ charges of the model multiplets. 



Multiplet 


X 





QlL 


QaL 




daR 


Ur 


DaR 




laR 


;B-charge 








1 
3 


1 

3 


1 
3 


1 
3 


1 

3 


1 
3 








£-charge 


4 
3 


2 

3 


2 
3 


2 
3 








-2 


2 


1 

3 


1 



where the subscripts p, m and n stand for SU(3)l indices. 

The VEV CO gives mass for the exotic quarks U, Da and the new gauge bosons Z' , X, Y, 
while the VEVs u and v give mass for the quarks Ua, da, the leptons la and all the ordinary 
gauge bosons Z, W [9]. To keep a consistency with the effective theory, the VEVs in this 
model have to satisfy the constraint 

■a^ < < uj^. (2.8) 

In addition we can derive v « 'Uwcak = 246 GeV and \u\ < 2.46 GeV from the mass of W 
boson and the p parameter [6], respectively. From atomic parity violation in cesium, the 
bound for the mass of new natural gauge boson is given by Mz' > 564 GeV (w > 1400 GeV) 
[10]. From the analysis on quark masses, higher values for oj can be required, for example, 
up to 10 TeV [9]. 

The Yukawa couplings of (2.6) possess an extra global symmetry [11, 12] which is not 
broken hy v,uj but by u. From these couplings, one can find the following lepton symmetry 
L as in Table 1 (only the fields with nonzero L are listed; all other fields have vanishing 
L). Here L is broken by u which is behind L{xi) = 2, i.e., u is a kind of the SLB scale. 

It is interesting that the exotic quarks also carry the lepton number (so-called lepto- 
quarks); therefore, this L obviously does not commute with the gauge symmetry. One 
can then construct a new conserved charge C through L by making a linear combination 
L = xT^ + yTs + CI. Applying L on a lepton triplet, the coefficients will be determined 

L = -Its + CI. (2.9) 

Another useful conserved charge B, which is exactly not broken by u, v and lo, is usual 
baryon number: B = BI. Both the charges C and B for the fermion and Higgs multiplets 
are listed in Table 2. 

Let us note that the Yukawa couplings of (2.7) conserve B, however, violate C with ±2 
units which implies that these interactions are much smaller than the first ones [9]: 
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In this model, the most general Higgs potential has very simple form [13] 

v{x, 4>) = Axh + iA4>H + Ai(x^x)' + M{(t>H? 

+A3(X^X)(0V) + Mix^m^x)- (2.11) 

It is noteworthy that V{x-,4') does not contain trilinear scalar couplings and conserves 
both the mentioned global symmetries, this makes the Higgs potential much simpler and 
discriminative from the previous ones of the 3-3-1 models [11, 12, 14]. The non-zero values 
of X &^cl (j) at the minimum value of V{x, (p) can be obtained by 



Xtx= T\"^^^2^ ^^^. (2.12) 



4A1A2 - A| " 2 
Xsul - 2Xinl _ 



4AiA2-Ai - 2- ^^-^^^ 



Any other choice of u, uj for the vacuum value of x satisfying (2.12) gives the same physics 
because it is related to (2.3) by an SU(3)l(X'U(1)x transformation. It is worth noting that 
the assumed u 7^ is therefore given in a general case. This model of course leads to the 
formation of Majoron [13]. 



3 Stable Higgs bosons in 3-3-1 model 

This section is to show that the economical 3 — 3 — 1 model furnishes a good candidate 
for self interacting dark matter. The important properties are that dark matter must be 
stable and neutral. Hence, we are going to consider the scalar sector of the model and 
specially neutral scalar sector, and we can specify whether any of them can satisfy the self 
interacting dark matter conditions. 

Let us review the Higgs states and coupling constants. In this model, the most general 
Higgs potential has very simple form given in (2.11). As usual, we first shift the Higgs 
fields as follows: 







^/2 1 




/ 




X = 


X2 












'72/ 









The subscript P denotes physical fields as in the usual treatment. The constraint equations 
at the tree level are given as 

^2 

/x? + Ai(n2 + a;2) + A3y = 0, (3.2) 
/.i + W + A3^^^^^ = 0. (3.3) 

Note that u is a parameter of lepton-number violation, therefore the terms linear in u 
violate the latter. Applying the constraint equations (3.2) and (3.3) we get the minimum 
value, mass terms, lepton-number conserving and violating interactions as follows 

V{X. 4>) = ^min + V^^ass + ^mass + ^LNC + ^InV, (3.4) 
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where 



A 



mm 
N 

mass 



2 4 
V 



1 



v 



C 



Ai(ti5i + wSa)^ + \2v'^Sl + A3?;(uS'i + <:j53)52, 

A4, 



-(U0J" +t;X2 +'^</'3 )Ml +^X2 +W'A3)> 



^LNC = \i{xh? + A2((/'V)' + A3(X^X)(</'V) + A4(xV)(0^x) 

+2Ai^^S3(x^x) + 2\2vS2{(t>^(t>) + A3t;S2(x^x) + A3u;53(0V) 



(3.5) 
(3.6) 



(3.7) 
(3.8) 



In the above equations, we have dropped the subscript P and used x = (xi5X2 ^Xs)"^) 



)-^. Moreover, we have expanded the neutral Higgs fields as 



Si + iAi 
Xi ^ 



^/2 ' 



53 + i^s 
X3 = 



,0 _ -5*2 + 

t>2 - 



V2 



(3.9) 



In the pseudoscalar sector, all the fields are Goldstone bosons: Gi = Ai, G2 = A2 and 
G3 = A^ (cl. Eq. (3.5)). The scalar fields S*!, S2 and 5*3 gain masses via (3.5), thus we get 
one Goldstone boson G4 and two neutral physical fields the standard model and the 
new H'^ with masses 



m]jo = \2V^ + Ai(n2 + _ ^[X2V^ _ x^{v? + a;2)]2 + Ai?;2(n2 + a;^) 

4A1A2 — Ao 2 

— ; -v\ 

2Ai 



(3.10) 



Mi 



'I, = X^v^ + Ai(n2 + + ^[A2t;2 _ Xiiu^ + a;2)]2 + A2^2(^2 + ^2) 



2Aiw2. 



In term of original fields, the Goldstone and Higgs fields are given by 

1 



G4 



1 + ti 



--{Si — teSs), 



i/o = c^S2 - -j^^iteSi + 53), 



= S(^S2 + 



1 + i 



-.{teSi + Ss), 



where 



t2C 



AiM2 - A2M2 ■ 



(3.11) 

(3.12) 
(3.13) 
(3.14) 



(3.15) 
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From Eq. (3.11), it follows that mass of the new Higgs boson M^jo is related to mass of 
the bilepton gauge (or via the law of Pythagoras) through 



X 



l + O 



Ml 



2\is 



w 



M 



X 



l + O 



Ml 



18.8AiAf 



Here, we have used a = and = 0.231. 



In the charged Higgs sector, the mass terms for 

ical seal 

m = ^^=^^(n(/.+ + vxt + ^4>t) 



(3.16) 

iiX2, (t>-i) are given by (3.6), thus 



there are two Goldstone bosons and one physical scalar field: 

1 



^/ V? -\- v"^ + Oj"^ 



with mass 



M: 



HJ 2 



— (U +V +UJ) 



2\. 



Ml 
9^ 



^ ^7\ )\^2 

27ra 



■Mp ~ 4.7A4M^. 



The two remaining Goldstone bosons are 

1 . ,^ 



3 !■> 



(l + t2)(n2+t;2+^2 



[t;(i,(/)+ + (/>+) - io{l + t^)x^; 



(3.17) 

(3.18) 

(3.19) 
(3.20) 



Thus, all the pseudoscalars are eigenstates and massless (Goldstone). Other fields are 
related to the scalars in the weak basis by the linear transformations: 



(3.21) 













h 




CqS0 






I 


ce 




cj^ + civ 



.2n,2 



UJS0 



vce 



ujce 



2 



-UJ 



VCq 





(4>t\ 


) 


xt 







(3.22) 



Let us comment on our physical Higgs bosons. In the effective approximation w ^ v,u, 
from Eqs. (3.21), and (3.22) it follows that 



5*2, Hi ~ 53, 



Ga 



Si, 

xt 



(3.23) 



From the Higgs gauge interactions given in [13], the coupling constants of Higgs and SM 
gauge bosons depend on sr with t2c = . In the w ^ v,u limit, Mx ^ Mw or 



l^2cl ~^ 0. Therefore, the Higgs does not interact with the SM gauge bosons , , 7. 
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However, there are couplings of Higgs with the Bilepton Y and Z' . In order to forbid 
the decay of Hf, we assume that M^o <M^. It means that 2XiUJ^ < jQ^uj'^ or Ai < 0.051. 
The interactions of Higgs with new gauge boson Z' is Z' — — G3 interaction. But G3 
is a Goldstone bosons, this interaction can be gauged away by a unitary transformation. 

Let us consider the interaction of dark matter to Higgs bosons. From the Higgs po- 
tential (2.11), we can obtain the coupling of the new Higgs H'l to H^H^ . The decay rate 




of the i?? is written as 



^m^H'^H^ = 16^ I ^ - 1 • ^^-24) 



The lifetime is the inversion of decay rate r = ^, with h = 6.6 x 10 GeV x s. If 
taking r > 10^*^5 (the life time longer than our universe's age), M^jo = 7000 GeV, tt; = 
10000 GeV, Mho = 120 GeV, then we get 



A^ lOW 2 X 120^ , ^ . 



rffO_^„or.o = ^ 1 ^ ~ 284 xAl (3.25) 

H-,^H H levrT.lO'^V 49.106 J 3 ^ ^ 

In order to get the constraint on the lifetime of longer than our universe's age, it is 
easy to see that the value of A3 is approximately order of 10~^^. It is to be emphasized 
that the limit of A3 makes sure that is small. 

To avoid decaying to H2H2, we need the constraint for the mass of two Higgs, 
namely M?,o < 4M^_,., which means Ai < A4. From the Lagrangian given in (2.5), it is 

easy to see that the Hf does not interact with the SM leptons but it interacts with exotic 
quarks. As we know the exotic quarks are heavy ones, we assume that their masses are 
heavier than that of Hf. Hence, Hf can be stable and be candidate for dark matter. 



4 Thermal relic abundance 
4.1 Constraints 

Before considering the relic abundance of dark matter, let us summarize the constraints 
on the couplings Ai^2,3,45 the VEV w, and exotic quarks masses: 

1. From Eqs. (3.10) and (3.11), we obtain the constraints as follows 

Ai > 0, A2 > 0, 4A1A2 > Ai. (4.1) 

2. The mass of the charged Higgs boson Hf' is proportional to that of the charged 
bilepton Y through a coefficient of Higgs self-interaction A4 > 0. Analogously, this 
happens for the standard-model-like Higgs boson and the new Hf . Combining 
(4.1) with the constraint equations (3.2), (3.3) we get a consequence: A3 is negative 
(A3 < 0). 
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3. In order to get the stable Higgs particle H^, we need the constraints as follows 

Ai < A4, Ai < 0.051, IA3I ~ 10-24, M^o < Mu. (4.2) 
Since A3 < 0, we get A3 ~ —10-^^. 

4. In the limit of A given in (4.2), the SM Higgs mass can be estimated as M'^q = 

—^^Xr~^^'^ ~ 2A2i'2. Combining with the constraint 80 < M^jo <160 GeV, we can 
obtain the constraints on A2 as follows: 0.053 < A2 < 0.212. 

4.2 Implication for parameter space from the WMAP constraints. 

In this subsection, we discuss constraints on the parameter space of the 3-3-1 model origi- 
nating from the WMAP results on dark matter relic density [15], 

nh"^ = 0.1120 ±0.0056 . 

In order to calculate the relic density, we use micrOMEGAs 2.4 [16] after implementing new 
model files into CalcHEP [17]. The parameters of our model are the self-Higgs couplings, 
Ai,A2,A3,A4, the VEV w and exotic quarks masses. Note that the usual quarks u,s,b 
gain masses at one- loop level [9]. We express the couplings of Higgs with exotic quarks 
s",S22,'S33 as functions of Ai, w and exotic quarks masses. 

All Feynman diagrams contributing to the annihilation of Hf Higgs are listed in the 
Appendix (A). At the tree level, the annihilation of dark matter can be done through s-, 
or t-channel, or direct annihilation. Since there is neither coupling of H^Hf to one neural 
gauge boson nor coupling of HfH^ with one fermion, the propagator in the s-channel can 
be or Higgs only. To draw Feynman diagrams contributing to the annihilation of 
HfH^ through s- and t-channels, we list all non-zero couplings H^AB and H^AB, where 
A, B can be Higgs, or gauge boson, or fermion. We see that couples to one usual 
quark (anti-usual quark) u, s, b and one anti-exotic quark (exotic quark) while couples 
to c, d, t quarks. Therefore, the annihilation of HIH^ into uu, ss, hh are done through 
t-channel through exotic quark exchange while the contributions of the remaining usual 
quarks are done via s-channel through exchange. Since the coupling H^H^H^ ~ vX^, 
the contribution of cc, dd, tt channels to is very small. Higgs can also annihilate 
into two Higgs bosons or two gauge bosons directly. Theses vertices arise from the Higgs 
potential and Higgs- gauge interactions. 

First we study the behavior of Qh'^ as a function of one parameter each time. Table 
3 shows the dependence of Qh^ on A2, A3, A4 corresponding to the point 1, 2, 3. In all 
three cases, we fix Ai = 0.04, w = W TeV, and exotic quarks masses Mu = 36 TeV, 
Mz)2 = -^Da = 100 TeV as a special choice of parameters in the WMAP allowed band 
(please look at Fig. 1), the green dot-dashed line. We can see that neither Qh"^ nor 
the contribution of channels change when varying A2 in the range 0.053 ~ 0.212, or A3 
from —10"^^ to —W~^^, which regions satisfy the constraints given in 4.1. The couplings 
H^H^H'^ and H^H2 H2 are proportional to f A2, and these contribute to the annihilation 
of dark matter through s-channel H exchange. The coupling HfH^H^ is vX^ ~ A3, 
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Table 3. ilh'^ and dominant channels when varying A2,A3,A4. 



point 


1 


2 


3 


Ai 


0.04 


0.04 


0.04 


A2 


0.053 to 0.212 


0.12 


0.12 


A3 


_10-24 


-10-33 to -10-20 


_10-24 


A4 


0.06 


0.06 


0.004 to 200 


w (TeV) 


10 


10 


10 


Mu (TeV) 


36 


36 


36 




100 


100 


100 


M^aCTeV) 


100 


100 


100 




0.1127 


0.1127 


(0.1116, 0.1128) 




nn(97.40%) 


uu{97A0%) 


uu{97.36 - 97.45%) 




ss(1.26%) 


ss(1.26%) 


ss(1.26%) 


channels 


66(1.28%) 


66(1.28%) 


66(1.28%) 




H+H2 {0.05%) 


H^H^ (0.05%) 


H+H^ (0-0.09%) 




rest (0.01%) 


rest (0.01%) 


rest (0.01 - 0.02%) 



where A3 is very small. That is why Qh"^ does not depend on A2. The relic density changes 
negligibly if we vary A4. It gets the minimum value O/i^ = 0.1116 at the point A4 = 0.15 
and the maximum value Qh'^ = 0.1128 at the two points A4 = 0.08 and A4 = 0.19. With 
A4 > 0.2, the relic density keeps constant value Qh'^ = 0.1095. 

Now we consider the dependence of the relic density of dark matter on the remain- 
ing parameters Xi^w, Mu, M^^ and Mjj^. First, we fix the values of A2,3,4 satisfying the 
constraints given in (4.2), especially taking A2 = 0.12, A3 = — lO-^^, A4 = 0.06 and varying 
the masses of exotic quarks. We consider the relic density as a function of Ai. Fig. 1 com- 
pares WMAP data to the theoretical prediction. The red dashed line presents predictions 
from our theory by fixing = M^g = 100 TeV, w=10 TeV and Mu = 24 TeV. In order 
to meet fully the WMAP dada, the value of Ai must be different from the allowed value 
in (4.2). However, if we change the masses of exotic quarks, we can obtain allowed region, 
namely the green dot-dashed line given by taking w = 10 TeV and Mjj = 36 TeV. The 
allowed region of Ai satisfy both the WMAP data and the stable Higgs constraints (4.2) is 
0.0393 < Ai < 0.0406. The orange full line is obtained by fixing u; = 30 TeV and Mu = 36 
TeV. In this case, the constraints on Ai is 0.0424 < Ai < 0.0436. On the other hand, if 
we change the mass of exotic D-quarks, we can find the other allowed region of Ai. For 
example, if we take Mjj^ = Mo^ = 12 TeV, the allowed region of Ai is 0.0502 < Ai < 
0.051. Hence, we could conclude that the mass of exotic [/-quark can be larger or smaller 
than that of D-quarks in order to come to agreement with the WMAP data. 

Fig. 2 shows the dependence of the relic density on the VEV w for Ai = 0.04, A2 = 
0.12, A3 = -10-24 and A4 = 0.06. This figure shows that the VEV w < 15.33 TeV is in the 
WMAP-allowed region for Mu = 36 TeV, Md2= Md^ =100 TeV. However, if the values of 
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0.14 


1 1 1 1 1 1 


,' ' \ 


0.13 




_ 

=0.051 


0.12 


f!/r=0.1176 




n 1 1 

^0.11 




in . ^' N \ 


0.10 
0.09 


fl/f=0.1064 


/ / — V- — 

' ;V ^^'^ - 


0.08 


" 





0.00 0.02 0.04 0.06 0.08 0.10 

Figure 1. nh'^ vs Ai for A2 = 0.12, A3 = -10^24^ ^ q.OG, Md^ ^Mds^IQQ TcV, and for w^lO 
TeV, Mu = 24 TeV (red dashed line), w = 10 TeV, Mu = 36 TeV (green dot-dashed Hue), w = 
30 TeV, Mu = 36 TeV (orange full line), and Md^ = Md^ = 12 TeV, w = 10 TeV, Mu = 70 TeV 
(brown large dashing line). The blue dotted vertical line corresponds to Ai = 0.051. 



Mu = 24 TeV or Mq = Mjj = 36 TeV, there is no allowed region of uj in agreement with 
the WMAP data. The situation becomes totally different for Mu = 70 TeV, M^a = Md^ 
= 12 TeV (brown large dashing line). The relic density at first increases then decreases as 
a function of w. In the WMAP band, w is in the range 8.752 - 13.85 GeV or 23.3 - 24.61 
GeV. 

Similarly, we can figure out the region of Mu, M^^, and M^,^ in agreement with the 
WMAP data by fixing Ai = 0.04, A2 = 0.12, A3 = -10~24^ A4 = 0.06, and u; = 10 TeV. The 
value of Mu is in the narrow band 34.93 ~ 36.73 TeV for Af^^ = Md^ = 100 TeV, and Mu 
is should be heavier, 66.71 < Mu < 85.01 TeV if we take M^^ = M^g = 12 TeV. In case 
of Mu = 36 TeV and Md^ = 100 TeV, Md2 is in the region of 37.99 < Md^ < 259.7 TeV. 
For Mu = 36 TeV and M^,^ = 100 TeV, the relic density is always in the WMAP-allowed 
region for Mjj^ > 40 TeV. There are many choices of {Mu, Md^, Md^) set satisfying the 
WMAP result, and we can see that the WMAP constraints do not require the order of 
exotic quarks masses. 

To give an overview of the behavior of the relic density in the M^^ - M^^ plane, as 
shown in Fig. 3, we consider the model with Ai = 0.04, A2= 0.12, A3 = — 10~^^, A4 = 0.06 
and "u; = 10 TeV. The region of Md^ - Md.^ in agrement with the WMAP is very wide 
for Mu = 36 TeV (red), while it seems to be two narrow bands for Mu = 70 TeV (grey 
bands). 

From now on we take M^ = Mu,^ = Md^ for convenience. In the Mu - Md plane (see 
Fig. 4), we can see that to satisfy the WMAP constraints, Mu must be heavier than 35.2 
TeV and Md must be heavier than 11.8 TeV. For Mu = 36 TeV and u; = 10 TeV, the 
relic density as slowly varying function of Md . The situation is similar to that of Mu = 40 
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Figure 2. Q.h^ vs w for Ai = 0.04, A2 = 0.12, A3 = -10^24^ A4 = 0.06, and for Mu ^ 24 TcV, Md.^ 
= AfD3 = 100 TeV (red dashed line), Mu = 36 TcV, A/d, = = 100 TeV (green dot-dashed 

line), Mu = Md^ = Md, = 36 TeV (orange full line) and Mu = 70 TeV, Md^ = Md^ = 12 TeV 
(brown large dashing line). 
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Figure 3. Contour plots for 0.1064 < Qli^ < 0.1176 (WMAP constraints) in Md^ - Md^ plane for 
Ai = 0.04, A2 = 0.12, A3 = -10-24^ A4 0.06, w=10 TeV, and for Mu = 36 TeV (red) and Mu = 
70 TeV (grey). 



TeV and ^« = 30 TeV. On the other hand, for the value of M/5 around 12 TeV, the rehc 
density varies very slowly as a function of Mu. For 12 < Mo < 22 TeV, the relic density 
at If = 10 TeV is the same as that a.t w = 2>{) TeV. 
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Figure 4. Contour plots for 0.1064 < VlK^ < 0.1176 (WMAP constraints) in Mu - Md plane for 
Ai=0.04, A2=0.12, A3 = -10^24^ A4=0.06 and for w=10 TcV (red) and ■u;==30 TcV (green). 



Similarly we investigate contour plots for 0.1064 < Q.h? < 0.1176 (WMAP constraints) 
in the Mu - w plane. The allowed region of Mjj is in the narrow band for fixing Md. For 
an example, if 15 < A/d < 21 TeV, the allowed region of Mu is in 40 < Mu < 64 TeV. 
Finally we study contour plots for 0.1064 < flh'^ < 0.1176 (WMAP constraints) in the Ai - 
Mu plane. Combining with the constraint on Ai given in (4.2), we can see that the allowed 
bands are: 0.028 < Ai < 0.051 for = Md3=W0 TeV, and for Md^ = Mds=12 TeV, 
0.039 < Ai < 0.051 if u; = 10 TeV, and no region of Ai allowed if w = 30 TeV. 

In next section we will study how to search for the DM in the WMAP - allowed region, 
in direct and indirect searches. We would like to analyze the results as functions of M^jo, 
which is expressed in term of Ai and w. With the WMAP constrains, we use the best 
parameter space, Ai = 0.04, A2 = 0.12, A3 = -10-^4^ A4 = 0.06, A% = 36 TeV, Md^ = 
Mu).^ = 100 TeV, and we vary 5 < w < 15.3 TeV, which requires M^^o to be few TeVs. 

5 Direct and indirect searches for the dark matter 

Experimentalists worldwide are actively chasing searches for DM candidates either directly 
through detection of elastic scattering of the weakly interacting massive particles with the 
nuclei in a large detector or indirectly through detection of products of the dark matter 
annihilation (photons, positrons, neutrinos or antiprotons) in the galaxy or in the sun. 

5.1 Direct search 

In direct search, the recoil energy deposited by scattering of WIMPs with the nuclei is 
measured. In general, WIMP-nuclei interactions can be split into two types: a spin inde- 
pendent interaction and a spin dependent interaction. In our model, scalar Hf Higgs DM 
can only contribute to spin independent interaction. To calculate the direct detection rate 
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we use the method, as mentioned in [18]; The direct detection rate depends on the WIMP 
nucleus cross section. To derive the ffj'-nucleus cross section one has to compute first the 
interactions at the quark level then convert them into effective couplings of WIMPs to 
protons and neutrons. Finally, we have to sum the proton and neutron contributions and 
turn this summation into a cross section at the nuclear level. 

The calculation of the cross section for WIMP scattering on a nucleon is obtained at 
the tree level. The normalized cross section on a point-like nucleus is given as 



^f^lo(^Zfp + iA-Z)U) 



TT 



(5.1) 



where Hjjo is the - nucleus reduced mass, fp and /„, are amplitudes for protons and 
neutrons, respectively. For Xenon, A = 131, Z = 54, while for Germanium A = 76, Z = 
32. The recoil spectrum of the nuclei depends on the velocity distribution and is contained 
in the elastic form factor of the nucleus. Using micrOMEGAs 2.4, we get the amplitudes 
and cross sections for WIMP-nucleon elastic scattering calculated at zero momentum as 
well as the total number of events/day/kg if we consider detector made of Xe or Ge. 
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Figure 5. DM-nucleon cross sections vs Af^o for Ai = 0.04, A2 = 0.12, A3 
0.06, Mu ^ 36 TeV and Md^ = Md^ = 100 TeV.' 



= -10-24, A4 = 



Fig. 5 shows the values of o"wiMP-nucieon as a function of dark matter mass by fixing 
the nucleon form factors, (Tq = 40 MeV and a-^N = 55 MeV. The value of o-wiMP-nucicon 
is in order of W~^{pb), which is allowed by experimental constraints of CDMs 2009 (Ge). 
However, in the limit the dark matter mass is smaller than 2.5 TeV or larger than 3.5 
TeV, the result given by our theoretical prediction is somehow different from experiment 
of CDMs 2009 (Ge). We would like to emphasize that the form factors of nucleons can 
be reset by changing the pion-nucleon sigma term, ct^at = 55-73 MeV and from the SU(3) 
symmetry breaking effect, uo = 35 it 5 MeV [19], however, the final WIMP-nucleon cross 
section predicted by our model does not change much. 
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Figure 6. Number of events/day/kg vs M^o for Ai = 0.04, A2 = 0.12, A3 = -lO'^^, A4 = 0.06, 
Mu = 36 TcV and M^, = M^,, = 100 TeV. ' 



Let us deal with the number of the nucleon recoil in Ge and Xe detectors. The predicted 
number of nucleon recoil is given in Fig. 6. The Xe detector is more sensitive than Ge 
detector. In the limit, 2.5 < M^jo < 3.5 TeV, the theoretical predictions are 22.3 and 14.6 
nucleon recoils those are observed in the Xe and Ge detectors for 1 kg per year, respectively. 
It is worth mentioning that the number of the nucleon recoils exposure between 10 keV ~ 
50 keV approximately equals one half of total of that number. 

5.2 Indirect search 

DM annihilation in the galactic halo produces pairs of the SM particles that hadronize and 
decay into stable particles. These particles then evolve freely in the interstellar medium. 
The final states with 7, e"*" and p are particularly interesting as they are the subject of 
indirect searches. From Feynman diagrams, we can see that the annihilation of H^H^ into 
tt^ cc^ dd, II, uivi and ZZ are done through s-channel exchange. The annihilation of 
H^Hf into H^H^ is done through s-, t-channel H'l, exchange or quartic couplings. 
Since the couphngs HIH^H^, HfH^H^, and H^H^H^H^ are proportional to A3, the 
contributions of those channels are very small. With the choice of parameters Ai = 0.04, 
A2 = 0.12, A3 = -10-24^ A4 = 0.06, Mu = 36 TeV, M^a = = 100 TeV and 5 < 

w < 15.3 TeV, the dominant channel is uu. For example, in case of w = 10 TeV, the 
relative contribution in % are displayed as following: (97.40% : H^Hf — )• uu); (1.28% : 
HfH^ bb); (1.26% : H^Hf ss); (0.05% : H^H^ H+H^); and (0.01% : the rest). 
The total annihilation cross section times the relative velocity of incoming dark matter 
particles is shown in Fig. 7. With the allowed region of the dark matter mass satisfied the 
WMAP constraints, we find 2.15 x 10~^^ < av < 2.4 x 10~^^ cm^/s. This result is in the 
same order of that given in [20] , which is said that away from the Higgs resonance and the 
W threshold, < cjf > is essentially constant and equal to the so-called typical annihilation 
cross section, < av >~ 3. 10~^^ cm^/s. The AMS-2 experiment given in [21] predicted for 
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dark matter mass up to 600 GeV. With dark matter mass from 100 GeV to 600 GeV, av 
keeps constant value in order of 10~'^^cm^/s. We expect that our result for heavy dark 
matter can be covered by the future experiment of AMS-2. 
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Figure 7. The annihilation cross section times the relative velocity of incoming DM particles vs 
Mho for Ai = 0.04, A2 = 0.12, A3 = -lO^^^^ A4 = 0.06, A% = 36 TeV and Md., = Md-, = 100 
TeV. 
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Figure 8. Photon flux vs E for Ai = 0.04, A2 = 0.12, A3 = -10^24^ ^ 0.O6, w = 10 TeV, Mu 
= 36 TeV and Md^ = Md^ = 100 TeV. 



Let us proceed with the discussion of the photon flux and energy spectrum for dark 
matter with mass 2828.4 GeV. The spectrum for photon flux is predicted in Fig. 8. It is 
easy to see that the photon flux in the energy range from a few MeV to 10 GeV is much 
larger than that of higher energy ranges. These results can be understood as following: As 
previously mentioned, our model predicts that the annihilation of the dark matter and anti- 
dark matter to u-quark and u-quark is the dominating channel. Therefore, the dominating 
jet is the neutral pion jet, composed of pairs of u-quark and anti-quark in this case. The 
7-rays from particle annihilation processes have spectra bounded by the rest mass energy 
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of the annihilation particle. The 7-rays are dominated by pion decay at low energy from a 
few MeV to 10 GeV. The additional contribution to photon spectrum at higher energy due 
to other annihilation processes such as polarization of the gauge bosons final state, photon 
radiation, etc, which are predicted to be tiny. Antiproton flux and positron flux also can 
be calculated by MicrOMEGAs 2.4. The fluxes go down fast as functions of energy and 
their values are significant at low energy as the same as photon case. 

6 Conclusions 

We have shown that the economical 3-3-1 model provides a candidate for dark matter 
without any discrete symmetry; and it just requires some constraints on Higgs coupling 
constant. The scalar Higgs Hf is a good candidate for self interacting dark matter. To 
forbid the decay of Hf, we require that Ai < 0.051, Ai < A4, and IA3I ~ 10~^^. The 
constraint on the Higgs coupling A3 looks unnatural, which could be canceled by introducing 
a discrete symmetry S3. However, by introducing new symmetry the Higgs sector becomes 
more complicated. Therefore, we do not consider such scheme in this work. The parameter 
space has been studied in detail and the results satisfying the WMAP observation are 
summarized as following: 

• The relic density does not change much when varying A2 from 0.053 to 0.212, A3 
around -lO'^^ value, A4 from 0.004 to 200. 

• Ai should be around 0.04. 

• The region of w is narrow compared to those of Mu , M^^ and . 

• [/-quark mass can be smaller or larger than D-quarks masses. 

We have studied direct and indirect searches for dark matter. The dark matter-nucleon 
cross section is in agreement with CDMs 2009 (Ge). The total number of events observed 
in Xe and Ge detectors is quite small because of the heavy dark matter. We hope that these 
results can be covered in future by experiments. Dark matter annihilation is considered 
with special choice of parameters. In case Mjj < M02 = M^^, the dominant channel is 
uu, while the dominant channel is ss for M^jj = M/jg > Mjj, because the interactions of 
quarks with other particles depend much on exotic quarks masses. However, choosing U- 
quark mass smaller or larger than D- quark mass does not affect our results on cross section 
times relative velocity as well as photon flux. The value of av is in order of 10~^^ cm^/s in 
agreement with typical annihilation cross section. Photon flux is dominated at low energy 
below 10 GeV. 
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A APPENDIX: Feynman diagrams contributing to the annihilation 
i^f? dark matter 
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